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Among systems that display generic scale invariance, those whose asymptotic properties are 
anisotropic in space (strong anisotropy, SA) have received relatively less attention, especially in 
the context of kinetic roughening for two-dimensional surfaces. This is in contrast with their ex- 
perimental ubiquity, e.g., in the context of thin-film production by diverse techniques. Based on 
exact results for integrable (linear) cases, here we formulate a SA Ansatz that, albeit equivalent to 
existing ones borrowed from equilibrium critical phenomena, is more naturally adapted to the type 
of observables that are measured in experiments on the dynamics of thin films, such as one and 
two-dimensional height structure factors. We test our Ansatz on a paradigmatic nonlinear stochas- 
tic equation displaying strong anisotropy like the Hwa-Kardar equation [Phys. Rev. Lett. 62, 1813 
(1989)], which was initially proposed to describe the interface dynamics of running sand piles. A 
very important role to elucidate its SA properties is played by an accurate (Gaussian) approximation 
through a nonlocal linear equation that shares the same asymptotic properties. 

PACS numbers: 05.40.-a, 68.35.Ct, 64.60. Ht 



I. INTRODUCTION 

Generic Scale Invariance (GSI) broadly makes refer- 
ence to the appearance of (quasi-) long-range order in 
sizable regions of the phase diagram of a given system. 
It can occur both in quantum and in classical systems 
and in equilibrium and out of equilibrium [ljj. In the 
case of classical, far-from-equilibrium phenomena, it pro- 
vides a natural framework for, among other, conspicuous 
systems such as those exhibiting self-organized-criticality 
(SOC) Thus, it embodies the behavior that was 

naturally attached to systems such as model sand piles, 
which were thought to reach a critical point without 
the need to fine-tune a control parameter. Although 
the mechanism for criticality of many SOC systems has 
turned out to be more elaborate than this Q, the GSI 
idea continues to play an important reference for non- 
equilibrium systems. Other specific contexts where it is 
at play are driven diffusive systems 0] and surface kinetic 
roughening [1]. 

Indeed, in terms of their natural description through 
Langevin equations with time-dependent noise the 
two main dynamical classes of kinctically rough surfaces 
(see, e.g., Ref. JM and references therein) necessarily 
display GSI 0, U; namely systems with nonconserved 
dynamics whose evolution equation is invariant under ar- 
bitrary shifts, h — >• h + const, in the values of the sur- 
face height h and systems with conserved dynamics and 
non-conserved noise. We will henceforth focus on these 
two groups of systems. A very important physical dif- 
ference with usual SOC systems concerns the relevant 
time scales [!, Thus, while in SOC systems there is 
a wide time-scale separation between driving (infinitely 
slow) and system response (infinitely fast), this is not 
the case for kinetically rough surfaces, for which noise is 



usually taken as time dependent, akin to thermal noise 
II- 

In the kinetic roughening context, to date most experi- 
mental and theoretical work has been devoted to charac- 
terize and understand scaling behavior that is isotropic in 
space. This is partly due to the fact that even the most 
natural anisotropic generalizations of the paradigmatic 
Kardar-Parisi-Zhang (KPZ) and conserved KPZ equa- 
tions, considered in Refs. and [Hj], respectively, turn 
out to display isotropic behavior themselves. However, to 
some extent, this should strike us as unexpected. Con- 
sidering, for instance, the application of kinetic rough- 
ening to standard contexts, like epitaxial growth of thin 
films [l3|], or crystal growth from a melt [3], physical 
anisotropics are ubiquitous in, e.g., energetic barriers for 
surfaces to relax by surface diffusion, or in surface tension 
effects. Thus, it would be natural to expect the occur- 
rence of anisotropic scaling properties in the asymptotic 
states of (some of) such systems. 

For instance, for thin films grown by Molecular Beam 
Epitaxy (MBE) on surfaces that are vicinal to a high- 
symmetry orientation, the direction of the average sur- 
face tilt and the average orientation of the ensuing steps 
play very different physical roles, frequently leading to 
anisotropic behavior as for the ste p b unching and/or 
meandering instabilities in Si(001) \mut- Even un- 
der morphologically stable MBE growth conditions, spa- 
tial anisotropies may occur, as for growth of GaAs films 
Qjl> [HI- Still with nanoscopic systems, erosion, rather 
than growth, of thin films by ion-beam sputtering (IBS) 
induces space anisotropies related with the different roles 
played by the direction on the target that lies along the 
projection of the ion beam and the direction perpendicu- 
lar to it [2(| . Reaching the realm of macroscopic spa- 
tially extended systems, fracture of solids provides still 
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another instance for the occurrence of space anisotropics, 
in this case between the crack propagation and crack 
front directions [HHl- 

Attempts have been made to explore the general role of 
anisotropics in kinetic roughening (see, e.g., Refs. 
and references therein), and some observations and char- 
acterizations of anisotropic behavior are available (26|. 
However, for the physical case of two-dimensional sur- 
faces, general conditions are not known on the occur- 
rence of isotropic vs anisotropic behavior, and neither is 
a formulation (Ansatz) of the latter available, that can be 
readily applied to analyze experimental data on, say, sur- 
face dynamics of thin films. Indeed, anisotropic kinetic 
roughening is encoded into a scaling Ansatz [27} that 
originates in the study of critical dynamics of equilib- 
rium statistical- mechanical systems [28| . Such a formu- 
lation is quite powerful from the theoretical point of view 
(enabling analysis of scaling properties for arbitrary sub- 
strate dimension, etc.) but is not particularly natural for 
the characterization of actual two-dimensional surfaces. 
In this way, the experimental community does not have a 
clear guiding principle that allows us to identify correctly 
the behavior that should be expected for each observable 
under conditions for strong anisotropy, and comparison 
with theoretical models is thus hampered. Alternative 
formulations are available [26| but have been tested onto 
linear models only [29j . so their generality has not been 
checked. Nor has a systematic study been done on the 
relation between different observables (like height corre- 
lation functions or ID and 2D height structure factors). 
Therefore, a fully consistent experimental analysis and a 
quantitative comparison with theoretical models cannot 
be done unless such relationship is fully clarified. 

The aim of this paper is to put forward a formulation 
of the dynamic scaling Ansatz for strongly anisotropic 
kinetically rough surfaces that is directly adapted to the 
analysis of experimental and/or simulation data for the 
physical case of two-dimensional interfaces, while being 
founded on the behavior displayed by model systems. 
To this end, we consider interface equations that can 
either be solved or that can be mapped into solvable 
cases. In general, we expect that progress achieved can 
allow a better assessment of the experimental impact of 
strongly anisotropic kinetic roughening, and, thus, add to 
an improved understanding of space anisotropies in large 
classes of GSI systems, for the physical two-dimensional 
case. 

In Scc.[H]of the paper, we, first, formulate our Ansatz 
and explore its implications for the various observables 
that are usually employed, both in experiments and in 
numerical studies of discrete or continuum models. The 
motivation for the Ansatz that is based on exact results 
for linear equations is discussed in Sec. IIIII Actually, re- 
cent experimental validation has been achieved in the 
interpretation of experiments of surface erosion by ion- 
beam sputtering [30(. Then, in Sec. IIV1 we again verify 
the Ansatz against a nonlinear equation that has played 
a relevant role in the context of space anisotropies. As 



a paradigmatic system, we specifically consider the Hwa- 
Kardar (HK) equation that describes the dynamics of a 
running sand pile [tl [3l[ and is a particular instance of 
conserved interface dynamics with nonconserved noise. 
Through direct numerical simulations, we analyze in de- 
tail its anisotropic roughening behavior in the light of our 
scaling Ansatz. Some apparent inconsistencies (sublead- 
ing corrections) are solved thanks to the mapping of the 
HK equation to a nonlocal linear equation that shares 
the same asymptotic properties (as conjectured HH). A 
discussion of these results is provided in Sec. [V] 

II. ANISOTROPIC SCALING ANSATZ 

For a two-dimensional rough interface, a straightfor- 
ward generalization of the height-difference correlation 
function that is usually considered for isotropic systems 

G(r)H[Mr + r )-Mro)] 2 ), W 

is provided by analogous correlation functions along each 
one of the two space directions, namely, 

G x (x) = ([h(xo + x,y ) ~ h(x ,y )} 2 ), (2) 

G y {y) = ([h(x ,y + y ) - h(x ,y )} 2 ). (3) 

In Eqs. ([1]) through ([3|), h(r) is the surface height above 
point r = (x, y) on a substrate plane where i"o = (xo, yo) 
is an arbitrary position, and brackets denote averages 
over the noise distribution (e.g., independent realizations 
of the experiment). In general, it is the physical con- 
ditions that dictate the two appropriate directions that 
play the role of x and y here. Thus, under IBS at oblique 
incidence, the projection of the ion beam on the target 
plane fixes one of them. Likewise, in MBE of vicinal sur- 
faces, it is the direction of the average tilt with respect to 
the height symmetry orientation that is fixing one main 
direction, and so on. 

While the statement of kinetic roughening in an 
isotropic system implies the simple power-law behav- 
ior G(r) ~ r 2a , where r = |r| and a is the roughness 
exponent the expected behavior in the presence of 
anisotropies is, rather [3lj . 

G x (x)~x 2a *, G y {y)~y 2a y. (4) 

where, in principle, two different roughness exponents 
exist, a x and a y . We will speak of strong anisotropy 
(SA) if they take different values in the asymptotic state 
of the system, while this will be referred to as of weak 
anisotropy (WA) if they take the same values. Note in- 
cidentally that in the present discussion we will assume 
that the surface is at a steady state at which statistical 
properties are time invariant. Later, we consider the time 
evolution explicitly 

Naturally, real-space correlations are not the only ob- 
servables one can measure from, e.g., simulations or from 
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actual experimental data. In the context of surface ki- 
netic roughening, a very important quantity is the two- 
dimensional power spectral density (PSD) or structure 
factor of the surface height, 



5(k) = (|Mk)p 



(5) 



where h(k) is the space Fourier transform of h(r)—h, with 
h being the space average of the height. On the one hand, 
this function can be easily measured in many experimen- 
tal setups by using, e.g., X-ray diffraction techniques [26| . 
On the other hand, for theoretical modeling and simula- 
tion, it has many advantages over real-space correlation 
functions, specially in the presence of crossover behavior 
and/or anomalous scaling; see, e.g., Ref. [l(| and refer- 
ences therein. As will be justified in the next sections, we 
will hypothesize the following long-time, large-distance 
behavior of the PSD for an anisotropic surface: 



5(k) =S{k Xl ky) 



vk 



2a,, 



(6) 



where the exponents measured in Fourier space a Xty will 
be related somehow to the roughness exponents a x , y and 
v is a mere constant. Our aim is precisely to see how 
the behavior we have just assumed translates into scaling 
behavior for other observables in the system. Note that 
Eq. (J6j> is an analog for nonconserved noise fluctuations 
of the behavior found in driven-diffusive systems Q . 

As mentioned above, frequently, e.g., in mentioned 
MBE or IBS systems, physical properties and geomet- 
ric constraints dictate the appropriate choices for the x 
and y directions. Nevertheless, as will be shown below, 
under conditions for strong anisotropy any choice of two 
orthogonal directions will lead to the same set of two dif- 
ferent exponents a x ,y, which guarantees the generality 
of Ansatz ([5]). In the case of fracture experiments, al- 
ternative choices for anisotropic scaling Ansatze are also 
available, in which, e.g., either an auxiliary dynamics is 
postulated [23|, [32| or expansions of observables over ap- 
propriate functional bases are performed that exploit the 
fact that isotropic materials often have anisotropic frac- 
ture surfaces only because of the breaking of isotropy by 
the initial conditions [3^, [34| . 

Once a choice of appropriate or convenient x, y di- 
rections has been made, two rather natural observables 
are the PSDs S x (k x ) and S y (k y ) of corresponding one- 
dimensional profiles. In e.g. thin film experiments, such 
observables are readily measured using scanning probe 
microscopies (AFM or STM) [26|. Thus, for instance, 
considering a fixed value y = yo, one defines 



S%\kx) — {h"yo\k'x^)h'yo\ ^a;))i 



where h yo (k x ) is the Fourier transform of the correspond- 
ing one-dimensional profile h(x,yo), 



hy (k x ) — 



1 



dx e lxkns h(x,y Q ) 



(8) 



Here we have considered the L — > oo limit of a finite sys- 
tem with lateral size L and periodic boundary conditions. 
Such a limit will be reconsidered when needed below. 

It is straightforward to obtain the relation between 
these observables [the same procedure as applied to 
S x (k x ) leads to the definition of S y (k y ) after exchang- 
ing labels x <-> y] and the two-dimensional PSD. Using 
([8]) in ([7]). and substituting h(x,yo) with its representa- 
tion in terms of the two-dimensional Fourier transform 
h(k x ,ky), 



h yo (k x ) 



1 



'2-k J- 



+ 00 



dk y e- lVok yh{k x ,k y ), (9) 



we obtain 



S x {k x ) — — 
Zir 

1 

7T 



+00 



dk y dk' y e-^^+^ihiK.kyyhi-K.k'y)) 

dk y S{k Xl k y ^j. 



Following the same steps, we obtain an analogous expres- 
sion for S y (k y ); in summary, 



£>x{k x ) — / dky S(k x ,ky), 



(10) 



(11) 



Now, substituting the scaling hypothesis ((6]) into Eqs. 
(fTOj) and we obtain the following scaling laws for 

the Fourier transform of the one-dimensional cuts: 



S y {k y ) — / dk x S(k x , fcy). 



S x {k x ~) 

Sy(ky) 



-(2a, 



l--( 2 «* 



-0 
-i/O 



(12) 
(13) 



where C = c\x/&y is the so-called anisotropy exponent. 
SA holds whenever £ ^ 1, whereas £ = 1 denotes WA, 
meaning that the steady state of the system is actually 
isotropic. As shown below (3f|, Eqs. (fT2"]) and (fTB")) arc 
equivalent to 



S x (kx) ~ k x ^ x ^ \ 
Sy(ky) ~ ky ^ v+ \ 



(14) 
(15) 



that provide the natural generalization to the SA case 
of the scaling behavior of the PSD of ID cuts of the 
surface in the isotropic case, in which a x — a y = a and 



S 



i ■ - 1 k x jy ^ j36j] . 
We still have to prove that the anisotropy condi- 
tion ( ^ 1 indeed corresponds to a x ^ a y , for which 
(7) we need to relate these real-space exponents with their 



momentum-space counterparts, a x . y . Such a relation is 
obtained from the definition of G x>y , Eqs. © and ([3]). 
The Fourier representation of the one-dimensional cut is 



h(x,y ) 



1 



+00 



(k x )e 



(16) 
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where h yo (k x ) is given by Eq. <j8j) . Using the latter into 
@, it is easy to see that 

2 f°° 

G x (x) = - / dk x [1 - cos(k x x)] S x {k x ), (17) 



a ID version of the general 2D relation between the 
height-difference correlation function and the height 
structure factor [26|. Now, using the scaling Ansatz for 
the one-dimensional PSD, Eqs. (fT2"j) and (fT5|) . we are able 
to write the scaling behavior of 



G x (x) — X 



2q x -C-1 



and 



G v (y) ~ V 



2c6 y -l/C-l 



2q 
2q 



25v 



l/C-i. 



(18) 

(19) 

(20) 
(21) 



Equations ([20]) and (|21j) lead immediately to Eqs. (fT4j) 
and p5[) and to the fact that a single anisotropy exponent 
is defined, namely 



c 



Ox 

a,. 



Oh. 

o ,, 



(22) 



In order to incorporate the time dependence into our 
scaling Ansatz, we consider another important observ- 
able, namely the global surface roughness. In anisotropic 
systems, it is natural to measure the roughness of slices 
of the surface along the x and y axes, W x (t) and W y (t), 
respectively. These are simply obtained as the averages 
of Eqs. ([2]) and ([3]) over all possible values of xo and j/o, 
respectively The Family- Vicsek (FV) dynamic scaling 
Ansatz [U that applies to isotropic kinetic roughening is 
typically formulated in terms of the short- and long-time 
behavior for the surface roughness W 2 (t) = ((h— h) 2 ) = 
J 5(k) dk. Thus @, W ~ for t < i 1/z , while W ~ L a 
for t 3> t x / z , where z is an independent exponent, t x l z is 
proportional to the length-scale below which nontrivial 
correlations have built up among height values at differ- 
ent substrate positions, and /? = a/z. 

In the SA case, one might expect that a similar be- 
havior would hold for the roughness functions of inde- 
pendent ID cuts, W XyV (t), but with two different growth 
exponents, j3 x = a x /z x and f3 y = a y /z y . However, it is 
not difficult to see that these two growth exponents coin- 
cide for any value of the anisotropy exponent £. To prove 
this, we consider two different anisotropic scale transfor- 
mations of the coordinates and the surface height, one 
along the x coordinate, r\ = (#i, j/i, 4i, /ii), and another 
one along the y direction, f*2 = (a?2, 2/2, ^2, /12), with 



ri = {bx,b c y,b z H,b a *h), 



r-2 



(&x,by,~b z yt,b a vh), 



(23) 
(24) 



where b and b are arbitrary positive factors. In these 
expressions, we have incorporated the shape of the dif- 
ferent rescaled coordinates under which scale invariancc 
holds. Notice, in principle one has to allow for a "re- 
sponse" of the time coordinate to a rescaling in space 
that is anisotropic, namely z x 7^ z y . Now the hypothesis 
of scale invariance implies that the same statistical prop- 
erties hold after we impose any of the two rescalings, (f2"5)) 
or (|24| (37( | . In particular, we can equate the two forms 
of rescaling. For the x coordinate, we get 6 = 6^. We 
can now apply the same procedure to the y coordinate 
and, from b^y = by, we get £ = C _1 - Then, by equating 
the two different ways of rescaling the time variable, from 
b Zx = b Zy we have ( = z y /z x . Finally, from the definition 



of the growth exponents j3 Xyl 



,y/ z x,y, we obtain 



Jy Z X Oty 



CC = i, 



(25) 



indeed implying /3 X = j3 y = j3 for any value of £. Sum- 
marizing, there are only three independent critical expo- 
nents characterizing a time-dependent SA surface, e.g., 
a x , C, and z x , from which all other exponents described 
in this section can be obtained. In particular, the behav- 
ior of the roughness of ID line profiles is W Xl y ~ t" for 
t < t 1 /*"'*, while W x , y ~ i Qx " for t > t 1 /**** with the 
corresponding values of exponents in terms of the three 
independent ones. In the next sections we will consider a 
number of examples for which we will fully verify all the 
scaling relations just derived. 



III. LINEAR EQUATIONS 

In this section we demonstrate a direct way of con- 
structing strongly anisotropic systems, by mixing two 
solvable one-dimensional equations with different scaling 
exponents. Thus we are able to combine two distinct 
physical phenomena acting along the two directions of 
the same system. As we will see, the resulting models 
will display SA that will be assessed through the scaling 
Ansatz proposed in the previous section. This will be 
done comparing the exact solution with direct numerical 
simulations, as a prelude to the next section in which a 
non-linear equation will be studied. 

For simplicity, the equation of motion will be written 
down for the Fourier components of the surface height. 
Thus, we consider 

fcftCM) = - (v x \k x \ n + v y \k y \ m )h(k,t) +fj(k,t). (26) 

Here k is defined in the two-dimensional space [— 7r,7r] 2 
and the term 7y(k, t) is the Fourier transform of a stan- 
dard time-dependent Gaussian white noise with zero 
mean and constant variance 2D. The equation just writ- 
ten actually defines a two-parameter family of models. 
Each single case corresponds to a specific choice of the 
pair of positive real numbers (n, m). In principle, here we 
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only consider positive values for coefficients u x ^ y so the 
dispersion relation of Eq. (j2"rJ)l , namely 

<r^) = -(y a \k ai \ n + v y \k y \ m ), (27) 

is negative definite. In this way, Eq. (|26|) is morphologi- 
cally stable. For related morp hologically unstable nonlin- 
ear models, see Refs. [38M4G |. Note also that, for generic 
values of n, m, Eq. (|26p is nonlocal in space, since it is 
not possible to inverse-Fourier transform it to have a lo- 
cal equation for the surface height h(x, y, t) in real space. 
This is only possible if n, m are even integers. See below 
for an example. 

Equation (|26|) can be explicitly solved for a flat initial 
condition ft,(k, t = 0) = 0, 

h(k,i)= [ dsexp[a(k)(t- s)]%(s). (28) 
Jo 

Using this solution, we can compute analytically the two- 
dimensional PSD for any of these linear models, 

D / 2v(k)t_ 1 \ 

Due to the morphological stability condition, the ex- 
ponential term vanishes in the t — > oo limit and Eq. 
([29]) reduces trivially to the scaling Ansatz (j6j. For a 
finite-size system, it is sufficient to verify the condition 
t > t sat = Ttiax X) y{L 1 / Zx ' v ) in order for the asymptotic 
form S{k x ,k y ) to approximate reasonably the full time- 
dependent S(k 7 t). 

To anticipate the type of analysis that will be per- 
formed in the next section for a more complex (nonlin- 
ear) equation, it is instructive to consider here a spe- 
cific example. Namely we consider a case that combines 
evaporation-condensation along the x axis and surface 
diffusion along the y axis. The resulting equation is a 
mixture of a ID Edward- Wilkinson (EW) equation Q 
along the former direction and a ID Linear Molecular 
Beam Epitaxy (LMBE) equation Q along the latter, 
namely 

d t h = vd 2 x h-JCd i y h + r ] {x 1 y 1 t). (30) 

This equation is simply the particular case (n, m) = (2, 4) 
of Eq. (|26p with v x = v and v y = /C, which can indeed 
be written down in real space, in terms of local operators 
acting on h(x,y,t). We will henceforth call this the 2-4 
equation, to recall the order of derivatives acting along 
each direction. Naturally, it is a particular case of a more 
general linear model in which evaporation-condensation 
and surface diffusion are fully anisotropic, namely, 

d t h=v x dlh + v y dlh r s 

4 2 2 4 ( 31 ) 

- K xx d x h - K, X yd x d y h - K vv d y h + i]. 

Equations of this form appear in the context of thin film 
growth and erosion [2^] and have partly been studied 
along the lines that follow. 
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FIG. 1: (Color online) One-dimensional projections S(k,0) 
(black circles, left axis) and S(0, k) (red diamonds, right axis) 
of the two-dimensional PSD of surfaces generated by the 2-4 
equation, (|30[) . Numerical integrations were performed with 
parameters v = 1, K, = 100, D = 1, L = 128, Ax = 1, and 
At = 1/2. The saturated two-dimensional PSD was averaged 
over 150 different noise realizations. The solid red line is a 
guide for the eye with slope —2, while the dashed blue line 
has slope —4. All units are arbitrary. 

Equation (|30|) being linear, we can indeed obtain all 
its scaling exponents exactly. Thus, the asymptotic two- 
dimensional PSD of the 2-4 equation behaves simply as 

S(k x ,ky)~ (isk 2 x +lCk±y\ (32) 

so a x = 1, iky = 2 and the anisotropic exponent is 
C = 1/2; hence, the system displays strong anisotropy 
as expected. The one-dimensional PSD of cuts along the 
coordinate axes should scale according to Eqs. (IT21 and 
P^)l . hence, 

S x {k x ) ~ k x , Sy(k y ) >~ k y , (33) 

from which we obtain direction-dependent roughness ex- 
ponents a x = 1/4 and a y = 1/2. Note [29[ that these 
values differ from those of the EW and the LMBE uni- 
versality classes in 1 and 2 dimensions Q. 

As far as the roughness is concerned, as we have seen 
before the growth exponent [3 is the same for the two 
directions. To compute its value, in this case it suffices 
to rescale Eq. ([30]) according to the transformation r\ ; see 
Eq. (|2"3"|) . Given the existence of the self-affine solution 
[Eqs. (|27l) and (|28)) ]. we can impose scale invariance on 
the rescaled equation, 

d t h = vb 2 - z *d 2 x h-/Cb 4 ^ z *d y 'h + b( 2a *- z * +1+ ^ 2 V , (34) 

namely one can request coefficients to be b indepen- 
dent. Thus, substituting the already-known exponents 
a x = 1/4 and £ = 1/2 we find the dynamic exponents 
z x = 2 and z y = 4, so j3 — 1/8. 
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FIG. 2: (Color online) One-dimensional cuts S(k cos 8, k sin 6) 
of the two-dimensional PSD of surfaces generated by the 2- 
4 equation (|30p along oblique directions in k space, 9 = 
0, 7r/12, 7r/6, 7r/4, and 7r/2, represented as functions of k. Nu- 
merical integrations were performed as in Fig. [1] Note the 
expected scalings —2 and —4 are best achieved for 6 close to 
and 7r/2, respectively. All units are arbitrary. 




FIG. 3: (Color online) PSD of one-dimensional cuts S x (k) 
(blue squares, left axis) and S y (k) (red circles, right axis) of 
surfaces generated by the 2-4 equation (|3Up . The numerical 
integrations were performed with the same parameters as in 
Fig. [T] The solid black line is a guide for the eye with slope 
—3/2, while the dashed green line has slope —2 and the dot- 
dashed magenta line has slope —4. All units are arbitrary. 



In order to check these analytical results, we have 
performed numerical simulations of the 2-4 equation by 
means of a pseudospectral integration algorithm as de- 
scribed in Ref. fdil ] and references therein. First, we have 
checked the validity of the asymptotic scaling Ansatz 
({52")) for the two-dimensional PSD. In Fig. Q] we present 
two projections of S{k Xl k y ) along the k x and k y axes, 
together with the analytical solutions, with very good 
agreement in the asymptotic regime. Note that, in this 
and all remaining figures, simulation data will be pro- 
vided through symbols, while theoretically expected scal- 
ing (whether approximate or exact) will be provided as 
reference lines. 

We can, moreover, verify that the same two exponent 
values a x y are actually obtained if one performs any cut 
of the PSD along any direction in the two-dimensional k 
space, see Fig. the x and y directions being optimal 
choices with respect to a clear-cut scaling behavior. This 
guarantees that the strongly anisotropic behavior found 
is not affected by an arbitrary choice of the x, y direc- 
tions in an associated physical system. Several different 
choices are shown in the figure, all of which take the form 
S(k cos 9, ksin9) for different values of the polar angle 9. 
Thus, e.g., S(k, 0), or S(0, k) in Fig. Q] are just the spe- 
cial 9 — and 9 = ir/2 cases, respectively. For other 
choices of 9, the scaling behavior of the corresponding 
cut S(k cos 9, k sin 9) crosses over between the expected 
k~ 2 and k~ A behaviors, as seen in Fig. [2] Again, from 
this point of view, 9 = 0,7r/2 allow us to elucidate the 
two existing independent exponents in the most unam- 
biguous way minimizing finite-size effects but do not play 
a more fundamental role. 

Further, we have also checked that the one-dimensional 



PSDs of cuts of the surface along the x and y directions 
are consistent with Eq. In Fig. [3] we see that indeed 
S x (k) (left axis) and S y (k) (right axis) display scaling 
behavior as predicted by the anisotropic Ansatz. Notice 
that, for small distances at which dynamics has not yet 
crossed over to its large-scale behavior, S y (k y ) ~ fc^ 4 , as 
would correspond for a ID profile of a ID LMBE system 

Finally, from the full dynamic evolution of the equation 
we obtained the value of the growth exponent (3. We 
have checked that this exponent does not depend on the 
direction chosen for its estimation. In Fig. |4] we show 
the roughness of the surface measured along each space 
direction. As expected, W x (t) ~ W y (t) ~ i 1 / 8 for times 
before saturation. 

In summary, for the 2-4 equation, all the anisotropic 
scaling relations derived in the previous section are ver- 
ified both numerically and analytically, thus providing a 
simple example of a system with strong anisotropy. How- 
ever, it would be interesting to probe these properties in 
a nontrivial system in which nonlincarities play a promi- 
nent role. The aim of the next section is to perform such 
a study for the celebrated example of the Hwa-Kardar 
equation. 



IV. NONLINEAR SYSTEMS: RUNNING SAND 
PILES 

The Hwa-Kardar equation (HK) was originally pro- 
posed in the context of self-organized criticality [3l[ . The 
intent was to construct a continuum field representa- 
tion of the sand-pile model proposed by Bak, Tang, and 
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FIG. 4: (Color online) Time evolution of the surface rough- 
ness functions measured along the two spatial directions, 
W x (t) (green circles, left axis) and W y (t) (blue diamonds, 
right axis), obtained from numerical integrations of the 2-4 
equation (|30[) . For visualization purposes, the value of W y (t) 
has been artificially offset vertically. The two solid black lines 
have the same slopes, 1/8. Parameters for the simulations are 
as in Fig. [T] All units are arbitrary. 



Wiesenfeld (BTW) @, @] and to study its critical behav- 
ior within such a "hydrodynamic" formulation. 

In the BTW model, sand grains are added randomly 
on a square lattice. At each lattice site, a discrete func- 
tion h keeps tracks of the height of the local sand col- 
umn. Whenever the local height value crosses a thresh- 
old, the sand grains of the columns are distributed among 
the nearest neighbors, until h becomes smaller than the 
threshold value. If one of the neighboring sites in turn 
crosses the threshold due to the grains received from the 
first site, the phenomenon propagates and triggers an 
avalanche. When an avalanche reaches the border of the 
lattice, the exceeding grains fall out of the system. Many 
versions of the this model can be found in the literature 
0, Q . In the original BTW model, the addition of each 
grain took place only when the system was completely in 
a relaxed state (without avalanches). Thus the addition 
of the grains (driving) and their transport (relaxation) 
occur at widely separated time scales. In Ref. [3l|, the 
authors focused on the opposite limit of a running sand- 
pile model, in which sand grains are added according to 
an external clock that is independent of the state of the 
system. 

Since transport of sand grains occurs only in the direc- 
tion parallel to gravity, the resulting field equation is in- 
herently anisotropic. Moreover, due to the fact that the 
relaxation mechanism implies a global loss of potential 
energy and the grains are allowed to escape through the 
boundaries without conservation constraints, the running 
sand-pile model is dissipative and open. The evolution 
equation is constructed taking into account the proper- 
tics of the discrete model jjj [3l|: The steady state is 



assumed to be on average a flat surface, h(r, t) measur- 
ing the local deviation from it. Here ru and rj_ are the 
components of r parallel and perpendicular to the trans- 
port direction, respectively. Because of the existence of a 
preferred flow direction, the system is not invariant with 
respect to the reflection symmetry in or h. However, 
it is assumed to be rotationally invariant in the rj_ di- 
rection, translationally invariant in both directions, and 
also invariant under the combined reflections h —h 
and ni — > — rii 0, [3l|. Taking into account all the model 
symmetries, and including the leading order nonlinearity, 
one obtains 

dth = v\\8f\h + vjy a ± h- -d\\h 2 + ri- (35) 

The first linear terms model the relaxation of the height 
through diffusive transport, whereas the nonlinearity ac- 
counts for the lack of inversion symmetry and is related 
to the presence of an external driving. In the absence of 
noise, dynamics is conserved, namely, Eq. (|35|) takes the 
form of a continuity equation, 

d t h + V • J = 0, (36) 

where 

J = -j/iVi/j- ^||«9||/i-^//jt. (37) 

Here t denotes the transport direction along the paral- 
lel coordinate. Given that the noise term r\ mimics the 
random addition of sand particles from outside the sys- 
tem, it is nonconserved, leading to GSI properties for the 
HK equation [H . Here, we will interpret these properties 
within our anisotropic scaling Ansatz, specifically for the 
HK equation in two dimensions (d = 2), 

d t h = v x dlh + v v dlh-^d x h 2 +r), (38) 

in which the x coordinate is parallel to the transport 
direction and the y coordinate is perpendicular to it. Al- 
though Eq. (|38p is nonlinear and cannot be solved ana- 
lytically, it is a very special system for which the three 
scaling exponents characterizing its SA behavior are be- 
lieved to be exactly known, due to the occurrence of large 
number of symmetries [H [3l[ . Specifically, these lead to 

z x = 2C, (39) 
z x + a x = l, (40) 
z x =2a x + (+l. (41) 

Equation (|39p originates from the fact that the nonlin- 
earity A, acting only along the x direction, cannot renor- 
malize linear operators that, like the one with coefficient 
v yi act along the perpendicular direction. Equation (|40[) 
implements a vertex cancellation that is believed to hold 
to arbitrary order in a pcrturbative expansion due to a 
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FIG. 5: (Color online) Top views of the morphology for sur- 
faces generated by use of the two-dimensional HK equation 
(|38|l . Numerical simulations have been performed with pa- 
rameters v x = v y = 1, A = 3, D = 1, L = 256, Arc = 1, 
and At = 10 -2 . The bottom panel has been obtained after 
rotating the top panel by 90°. 

concomitant symmetry of the equation under a Galilean- 
type transformation parametrized by A. Finally, Eq. (|4"Tj) 
is a hyperscaling relation (for d = 2) 0] deriving from the 
fact that Eq. (|38p implements conserved dynamics with 
nonconserved noise. As mentioned, the three scaling laws 
Eqs. (|3"9")) - (f4"Tj) are believed to hold exactly for arbitrary 
substrate dimension; see Refs. [HHlJ. For each value of 
d, such a nonhomogencous linear system actually has as 
a unique solution the set of three independent exponents 
that fully characterize the ensuing SA behavior. In our 
2D case, the values are 

a x = 4 «. = §, C =|. (42) 

From the relation z x = a x /(3, we get /3 = —1/6. Since 
the growth exponent is negative, subleading corrections 
appear in the time evolution of the roughness that com- 
plicate the estimation of j3 from the study of, e.g., W x {t) 



and W y (t). However, the scaling relations derived in Sec. 
HI! allow us to to compute 

Oly = --, Zy = 2, (43) 

and the values of the momentum-space exponents of the 
two-dimensional PSD S(k x ,k y ), 

3 

a x = -, & y = 1. (44) 

These will be simpler to check in simulations. The fact 
that all the critical exponent values can be calculated an- 
alytically makes the HK equation a perfect case study in 
which to check the scaling Ansatz (|6]) as a paradigm of 
a nonlinear system displaying strong anisotropy. More- 
over, to the best of our knowledge, to date there is no 
available study of the 2D HK equation through direct 
numerical simulations, in which the anisotropic scaling is 
analyzed. 

As for the 2-4 equation, the scaling behavior of the 2D 
HK equation (|38|) has been verified through pseudospec- 
tral numerical simulations. An example of the surface 
morphology that is obtained for long simulation times is 
provided in Fig. [5] In order to get a visual impression 
on the anisotropy of the system, in the same figure we 
present the same topography after a 90° rotation in the 
(x, y) plane. More quantitatively, in Fig.|5]wc present the 
two projections of the two-dimensional PSD S(k, 0) (left 
axis) and S(0, k) (right axis). Again, we clearly observe a 
scaling behavior like the one predicted by our anisotropic 
Ansatz, using the exact values Eq. (|4"4"|). 

However, for the PSDs of ID cuts of the surface, we 
find a striking difference. As can be seen in Fig. [3 while 
the long-distance behavior of S x {k x ) is as expected, this 
is not the case for the PSD of the one-dimensional cuts 
along the y direction, S y {k), that shows an unexpected 
behavior in the region of small wave vectors. The dis- 
crepancy between the numerical data and our theoretical 
prediction consists of two different effects: (i) The value 
of the PSD for the smallest value of k, i.e., S v (2tt/L), 
behaves quite differently with respect to the other points 
in the curve, "spoiling" the scaling behavior. Further 
numerical inspection (not shown) demonstrates that this 
feature is not due to fluctuations or numerical artifacts. 
(ii) Even neglecting this point, the scaling behavior of 
S y (k) does not resemble in any way the theoretical one 
given by Eqs. (|4"3")l . This unexpected behavior of S y (k) 
deserves to be studied with more detail since it could 
have important consequences about the anisotropic scal- 
ing Ansatz ([5]). As potential explanations, we can con- 
ceive of three possibilities: The Ansatz needs modifica- 
tions in the presence of nonlincaritics, the particular val- 
ues of the exponents induce this behavior, or both the 
nonlincarity and the exponent values concur in producing 
it. In order to elucidate this phenomenon, in the next sec- 
tion we construct a linear equation with the same scaling 
exponents of the HK equation. This nonlocal Gaussian 
approximation can be solved analytically and satisfies the 
anisotropic scaling Ansatz. 
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FIG. 6: (Color online) One-dimensional projections S(k, 0) 
(black circles, left axis) and S(0, k) (red diamonds, right axis) 
of the two-dimensional PSD for surfaces generated by the two- 
dimensional HK equation (|38l) . Numerical simulations have 
been performed with parameters as in Fig. [5] The saturated 
two-dimensional PSD was averaged over 100 different realiza- 
tions. The solid red line is a guide for the eye with slope 
—6/5, while the dashed blue line has slope —2. All units are 
arbitrary. 
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FIG. 7: (Color online) PSD of one-dimensional cuts S x (k) 
(blue squares, left axis) and S y (k) (red circles, right axis) 
of surfaces generated by the two-dimensional HK equation 
(|38[) . Numerical integrations were performed with the same 
parameters as in Fig. [6] The solid black line is a guide for 
the eye with slope — 2/3 while the dashed green line has slope 
— 1/3. Note that S y (k) does not follow the predicted scaling 
behavior at small k. All units are arbitrary. 



A. (Non-local) Gaussian approximation 

The fact that the HK asymptotics is, in principle, de- 
termined by three scaling relations allows us to put for- 
ward a linear equation that shares the same three inde- 
pendent exponents exactly. The key observation appears 
already in Eq. ([34)) of the previous section. Namely, 
for an (n, m) equation, a rescaling allows us to fix all 
three exponents by requesting statistical invariance. In 
order to see it, let us rescale Eq. ([26]) for n = 2a x , 
m = 2a yi and d = 2, according to the transformation 
hi = b a *h, k' x = b~ x k x , k' y = b'^ky, and t' = b z H 
[which is the momentum-space counterpart of r\ , see Eq. 
(|23j) ]. Requesting scale invariance, we obtain three rela- 
tions among the critical exponents, namely 

2a x = z x , (45) 
2(a y = z x , (46) 
z x = 2a x + C + 1. (47) 

From Eqs. ([20]) , ([21]) , and ([25]) , it is straightforward to see 
that the three relations (|45[) -(|47 ]) are completely equiva- 
lent, so only one of them is independent. Actually, Eq. 
(|47[) is hyperscaling, which holds for any linear equation 
and is also shared by the HK equation. For a noncon- 
served nonlinearity (for example, of the KPZ type), the 
noise term becomes renormalized and a such constraint 
may not hold. In such a case, the nonlinear equation and 
its Gaussian approximation would not share the same set 
of scaling exponents, in contrast with the case considered 
here. 



After these preliminary considerations, it is straight- 
forward to find a (n, m) equation that has the same 
three independent exponents as the nonlinear HK equa- 
tion. Thus, we just need to take n = 2a x = 6/5 and 
m = 2a y = 2. That is, the equation 

d t h(k,t) = - (is x \k x \ 6 / 5 + is y \k y \ 2 ) fc(k,t) + »KM), 

(48) 

has, by construction, the exponents given in ([42]) . Note 
that this equation is non-local in space since n = 6/5 
is not an even integer. Moreover, it is important to 
stress that Eq. (]48[) is not the linearized HK equa- 
tion. The latter would be obtained by setting A equal 
to zero, and, therefore, both exponents in momentum 
space would be 2a x = 2a y = 2, corresponding to 
the simple two-dimensional Edwards- Wilkinson equa- 
tion On the other hand, Eq. ([48]) can be consid- 
ered a Gaussian approximation of the HK equation, in 
the sense that it is a variational equation Q that has 
the asymptotic Boltzmann height distribution V{h} cx 
cxp[-(l/2L>)/(^|fc :c | 6 / 5 + ^|fc. y | 2 ) |/i(k)| 2 dk]. The fact 
that accurate approximations of nonlinear stochastic 
equations can be achieved through appropriate Gaus- 
sian counterparts is already known. Thus, in equilib- 
rium, the Bogoliubov-Feynman's inequality leads to a 
variational mean-field approximation that can be applied 
even to systems for which nonlinearities are non poly- 
nomial in the height field [42], [43|, while, out of equilib- 
rium, the so-called self-consistent expansion (see Ref. [44| 
for an overview with references) allows us to account for 
the scaling exponents of a number of-nontrivial systems 
through an appropriate Gaussian Ansatz. 
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FIG. 8: (Color online) One-dimensional projections S(k, 0) 
(black circles, left axis) and 5(0, k) (red diamonds, right 
axis) of the two-dimensional PSD of surfaces generated by 
the Gaussian approximation of the HK equation, Eq. (|48[1 . 
Numerical integrations have been performed with parameters 
v x = v y = 1, D = 1, L = 256, Ax = 1, and At = 0.1. The 
saturated two-dimensional PSD has been averaged over 300 
different realizations. The solid red line is a guide for the eye 
with slope —6/5 while the dashed blue line has slope —2. All 
units are arbitrary. 



Since all the scaling relations that we have derived stem 
from the two-dimensional PSD, for Eq. (fl"5)) we expect 
the behavior of S x (k) and S y (k) to be exactly as pre- 
dicted by CE) and (|T3]l. Before analyzing the PSD of 
one-dimensional cuts of the surface, we have checked the 
behavior of the two-dimensional PSD. In Fig. [5] we plot 
the projections of the 2D PSD onto the k x and k y axes, as 
in previous cases. The scaling for small k is as expected. 



Moreover, S x (k), reported in Fig. [9] (left axis), behaves 
as expected. However, the one-dimensional PSD of cuts 
along the y axis S y {k) (right axis) displays the same un- 
expected behavior as for the nonlinear HK equation, the 
first point being "shifted" above the rest of them and 
not adjusting to a straight line with the expected slope. 
Thus, such behavior does not seem due to the scaling 
Ansatz we are employing since, for a linear equation like 
(|48| , this simply expresses the exact behavior of the ob- 
servables in the long time limit. Thus, the behavior of 
S v {k) seems, rather, induced by the present exponent 
values. 



In order to clarify this issue, it is important to keep in 
mind that we are working for a finite system of lateral 
size L < oo, and with a finite lattice spacing a = Ax ^ 0. 



FIG. 9: (Color online) PSD of one-dimensional cuts S x (k) 
(blue squares, left axis) and S y (k) (red circles, right axis) of 
surfaces generated by the Gaussian approximation of the HK 
equation [Eq. (|48[l ]. The numerical integrations have been 
performed with the same parameters as in Fig. [8] The solid 
black line is a guide for the eye with slope —2/3 while the 
dashed green line has slope —1/3. Note that S y (k) does not 
follow the predicted scaling behavior. All units are arbitrary. 



The discrete counterparts of Eqs. (fT0|) and (JTT|) arc 
1 



S x \k x ) — 



Sy(ky) 



L 

1 

L 



N/2 

S(k x ,0) + 2 S(k x ,k v ) 

N/2 

S(0,ky)+2 ^ S(k X ,ky) 



, (49) 



■ (50) 



where = L/Ax = L/a, and n x ,n y = 1, . . . , N are dis- 
crete indices so k x ,y = 2Tm x ^ y / L. The hat stresses the 
discrete nature of these formulas. In the thermodynamic 
limit (for L — ¥ oo), these two series can be well approxi- 
mated as 



S x (k'x) ~ 



7r j a 



dky S{k x , ky) 



= k-^ a ^I(ki,a x ), 
1 f*/ a 

Sy{k y ) ~ / dk x S{k xi ky) 
n Jo 



where we have defined 



I(k,a) = - / — — ^. 

7T J 1 + u Za 



(51) 



(52) 



(53) 



In case the integral / is finite in the continuum a — > 
limit, the approximations (|5ip and (1521) hold and we in- 
deed obtain the scaling laws (fT2"|) and (|13|) . However, 
depending on the value of the critical exponents and the 
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FIG. 10: (Color online) Discrete ID PSDs S x {k) (left axis) 
and S y (k) (right axis) calculated from the sampled two- 
dimensional PSD S(k x ,k y ) of the Gaussian approximation 
Eq. 1)48 [I for increasing system size L. Other model parame- 
ters are the same as in Fig. [5] All units are arbitrary. 



FIG. 11: (Color online) Comparison of simulation data for 
S v (k) as obtained for Eq. ((48)) in Fig. [I0l with the analyt- 
ical solution of Eq. (|50[1 given by Eq. (|54p (solid blue line). 
The dashed red line corresponds to the asymptotic behavior 
S y (k) ~ fc -1,/3 . All units are arbitrary. 



numerical parameters employed, convergence of the inte- 
grals in (|5T|) and (|52)) can be too slow, and a clear-cut 
scaling may be jeopardized. 

In order to verify this issue, we have evaluated the 
series S x (k) and S y (k) from the known two-dimensional 
PSD S(k x , k v ) for increasing system size L (at fixed space 
resolution a). In Fig.[TU]we see that the asymptotic scal- 
ing is already clear-cut for S x {k) for the L values con- 
sidered, while for S y (k) it is still far from being seen. 
However, in the figure one can see that the curves for 
S y (k) become smoother as L increases, supporting the 
prediction that, for a sufficiently large system size, a well- 
defined scaling behavior could be seen. A further step 
can be taken by solving the integral in Eq. (f52|) and ([53)) 
exactly, 

lim S y (k y ) k~ 2& y a' 1 2 F 1 (p 1 ,l;l+p 1 ,p 2 ) , (54) 

where 2-P1 is a hypergeometric function with parameters 
pi = l/2a x and p 2 = -{tt / af^k' 2 ^ . In Fig. [UJ we 
show this approximation of S y (k y ) for a huge system 
(L = 10 8 ) with lattice cutoff a = 1. For the largest 
system sizes, the asymptotic scaling is clear and agrees 
perfectly with the one predicted from our scaling Ansatz. 
In the figure we, moreover, compare the analytical behav- 
ior with the numerical data from simulations of different 
system sizes. The numerical data indeed tend to adjust to 
the theoretical curve. This proves that, for the present 
values of the exponents, the sum ([50]) for S y (k y ) con- 
verges too slowly to the integral (|54j) for any numerically 
feasible system size. On the other hand, S x (k x ) displays 
the asymptotic behavior already at relatively small val- 
ues of L. Thus, the convergence rate strongly depends 
on the values of the critical exponents. One way to di- 
rectly verify this would be to perform numerical simu- 



lations with an adequate parameter set, corresponding 
to a theoretical curve in which the scaling behavior is 
clear. Unfortunately, as suggested by Fig. [TTJ this would 
require a computational effort beyond our capabilities. 

These results for the Gaussian approximation of the 
HK equation support the hypothesis of a similar slow- 
convergence problem also for the nonlinear HK equation 
proper. The fact that the behavior of S y (k y ) for Eqs. ([38)) 
and (|48p is qualitatively the same, and the fact that the 
presence of a nonlinearity usually contributes to slowing 
the convergence rate of observables, lead us to think that 
the asymptotic state for HK equation has not been fully 
reached in our simulations. However, as for its Gaussian 
approximation, we have not been able to perform a suffi- 
ciently large-scale numerical simulation due to the great 
computational effort required. 

The present slow-convergence problem may have im- 
portant implications for the interpretation of experimen- 
tal data. As pointed out earlier, many experimental se- 
tups are designed to measure the PSD of one-dimensional 
cuts, but we have shown here that for such integrated 
quantities there are systems in which the scaling may 
not be clearly detected. On the other hand, the two- 
dimensional PSD appears to be a more robust quantity, 
since it scales correctly already for relatively small val- 
ues of the system size L and relatively large values of the 
spatial discretization. We have recently confirmed these 
expectations in the analysis of experiments on surface 
erosion by ion-beam sputtering [30I ] . 



V. DISCUSSION AND CONCLUSIONS 

In this work we have presented a dynamic scaling 
Ansatz for two-dimensional anisotropic kinetic roughen- 
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ing that has a form readily applicable to the analysis 
of simulation or experimental data. The Ansatz incorpo- 
rates the behavior that characterizes strongly anisotropic 
linear equations, as has been thoroughly checked via their 
exact solution and numerical simulations. In this pro- 
cess, we have introduced a family of linear models that 
display SA while generally being nonlocal in space. Inci- 
dentally, nonlocal models (incorporating, e.g., fractional 
Laplacians) are lately finding ubiquitous use in equilib- 
rium and non-equilibrium statistical mechanics [451 ]; see 
Refs. [38|-|4(| and references therein for the specific con- 
text of kinetic roughening. 

A natural second step that we have taken in the present 
paper has been the validation of our scaling Ansatz 
against a numerical study of a paradigmatic nonlinear 
model of SA and in general of GSI systems, namely the 
HK equation for running sand piles. Again, agreement is 
quite satisfactory (as occurs for some experimental sys- 
tems [3(| ) and can be concluded after elucidation of slow 
convergence properties of certain observables, induced by 
(small) critical exponent values. This has required us 
to formulate and study a linear equation that has the 
same exact exponents as the HK nonlinear system and, 
thus, can be considered a Gaussian approximation of the 
latter. Thus, our results for a representative nontrivial 
anisotropic system suggest that the scarcity in assess- 
ments of anisotropic scaling properties in practical cases 
may be related with the occurrence of similar finite size- 
effects. 

Methodologically, our present work provides some re- 
sults that can be of interest in the analysis of strongly 
anisotropic interfaces. On the one hand, we have seen 
that in the case of systems with conserved dynamics and 
nonconserved noise it is always possible to tailor a fully 
dynamical description through an appropriate (n, m) lin- 
ear equation. This is due to the fact that the ensuing 
hypcrscaling relation reduces the number of independent 
exponents to two that are then readily tuned by choosing 
the real (non-necessarily integer) parameters n, m appro- 
priately. Such a procedure to put forward an effective dy- 
namical model may be of interest if dealing with complex 
nonlinear systems and / or with the description of simula- 
tion and/or experimental data. Incidentally, for noncon- 
served dynamics, some examples are also known on the 



equivalent stationary-state properties of some nonlinear 
models — such as the anisotropic KPZ equation for non- 
linearities with opposite signs — and of linear equations 
[46| . Likewise, a whole family on nonlocal and nonlinear 
models exists which shares the same probability distribu- 
tion function at steady state with nonlocal linear models 
akin to the present (n,m) equations; see Ref. [47| . 

On the other hand, our analysis has shown the prac- 
tical convenience of the use of certain observables (e.g., 
2D PSD function) over others (e.g., PSD functions of ID 
cuts of the 2D surface), provided one can access data 
with a large-enough signal-to-noise ratio. This has actu- 
ally been recently assessed in the analysis of experimental 
data from thin-film production [30| . 

In order to go beyond the present paper, one possibil- 
ity would be to consider the applicability of the exponent 
inequality that has been recently found for large classes 
of growth models in Ref. 48j, which leads to estimates 
of the scaling exponents even in nonlinear systems. This 
would require a generalization of the behavior of the var- 
ious observables involved (response and correlation func- 
tions) in the anisotropic case. Another natural avenue 
to push forward our present results is to consider the 
pertinence of our scaling Ansatz to nonlinear equations 
whose dynamics lies beyond reach of linear (Gaussian) 
approximations. Questions of interest in this connection 
are the occurrence, if at all, of SA behavior of the type 
we have discussed here and under which conditions. Note 
in particular, that all the systems we have analyzed so 
far correspond to conserved dynamics. Thus, the case of 
nonconserved system becomes particularly intriguing, es- 
pecially in view of the fact that the natural paradigmatic 
system, namely the anisotropic KPZ equation, does not 



show strongly anisotropic scaling in two dimensions [11 



Acknowledgments 

We are indebted to A. Keller for discussions. Partial 
support for this work was provided by MICINN (Spain) 
Grant No. FIS2009-12964-C05-01. E. V. acknowledges 
support by Universidad Carlos III dc Madrid. 



[1] D. Belitz, T. R. Kirkpatrick, and T. Vojta, Rev. Mod. 
Phys. 77, 579 (2005). 

[2] G. Grinstein, D.-H. Lee, and S. Sachdev, Phys. Rev. Lett. 
64, 1927 (1990). 

[3] G. Grinstein, J. Appl. Phys. 69, 5441 (1991). 

[4] G. Grinstein, in Scale Invariance, Interfaces, and Non- 
Equilibrium Dynamics, edited by A. J. McKane, M. Droz, 
J. Vannimenus, and D. Wolf (Plenum Press, New York, 
1995). 

[5] K. Christensen and N. R. Moloney, Complexity and Crit- 
icality (Imperial College Press, London, 2005). 



[6] R. Dickman, M. A. Munoz, A. Vespignani, and S. Zap- 
peri, Braz. J. Phys. 30, 27 (2000). 

[7] B. Schmittmann and R. K. P. Zia, Statistical mechan- 
ics of driven diffusive systems, in Phase Transitions and 
Critical Phenomena, Vol. 17, edited by C. Domb and J. 
L. Lebowitz (Academic Press, London, 2000). 

[8] A.-L. Barabasi and H. E. Stanley, Fractal Concepts 
in Surface Growth (Cambridge University Press, Cam- 
bridge, UK 1995). 

[9] M. Kardar, Statistical Physics of Fields (Cambridge Uni- 
versity Press, Cambridge, UK, 2007). 



13 



[10] R. Cuerno and L. Vazquez, in Advances in Statistical and 

Condensed Matter Physics, edited by E. Korutcheva and [29 

R. Cuerno, (Nova Science, New York, 2004). 
[11] D. E. Wolf, Phys. Rev. Lett. 67, 1783 (1991). [30 
[12] H. Kallabis, J. Phys. A: Math. Gen. 31, L581 (1998). 
[13] C. Misbah, O. Pierre-Louis, and Y. Saito, Rev. Mod. [31 

Phys. 82, 981 (2010). 
[14] S. H. Davis, Theory of Solidification (Cambridge Univer- [32 

sity Press, Cambridge, UK, 2001). 
[15] C. Schelling, G. Springholz, and F. Schaffler, Phys. Rev. [33 

Lett. 83, 995 (1999). 
[16] M. Sato and M. Uwaha, Phys. Rev. E 60, 7120 (1999). [34 
[17] A. D. Verga, Phys. Rev. B 80, 174115 (2009); 

arXiv.org: 1207.4354vl (2012). [35 
[18] A. Ballestad, B. J. Ruck, M. Adamcyk, T. Pinnington, 

and T. Tiedje, Phys. Rev. Lett. 86, 2377 (2001). [36 
[19] A. Ballestad, B. J. Ruck, J. H. Schmid, M. Adamcyk, 

E. Nodwell, C. Nicoll, and T. Tiedje, Phys. Rev. B 65, [37 

205302 (2002). 

[20] W. L. Chan and E. Chason, J. Appl. Phys. 101, 121301 [38 
(2007). 

[21] J. Munoz-Garcfa, L. Vazquez, R. Cuerno, J. A. Sanchez- [39 
Garcia, M. Castro, and R. Gago, in Toward Functional 
Nanomaterials, edited by Z. M. Wang (Springer, New [40 
York, 2009). 

[22] M. J. Alava, P. K. V. V. Nukala, and S. Zapperi, Adv. [41 
Phys. 55, 349 (2006). 

[23] D. Bonamy and E. Bouchaud, Phys. Rep. 498, 1 (2011). [42 

[24] U. C. Tauber and E. Frey, Europhys. Lett. 59, 655 

(2002). [43 

[25] U. C. Tauber, Critical Dynamics: A Field [44 
Theory Approach to Equilibrium and Non- 
Equilibrium Scaling Behavior, unpublished. [45 
http : //www. phys . vt . edu/-~tauber 

[26] Y. Zhao, G.-C. Wang, and T.-M. Lu, Characterization of 

Amorphous and Crystalline Rough Surface: Principles [46 
and Applications (Academic Press, San Diego, 2001). 

[27] B. Schmittmann, G. Pruessner, and H.-K. Janssen, Phys. [47 
Rev. E 73, 051603 (2006). [ 

[28] M. Henkel, H. Hinrichsen, and S. Liibeck, Non- 
Equilibrium Phase Transitions, Vol. I: Absorbing Phase 



Transitions (Springer, Dordrecht, 2008). 

Y.-P. Zhao, G.-C. Wang, and T.-M. Lu, Phys. Rev. B 

58, 13909 (1998). 

E. Vivo, M. Nicoli, M. Engler, T. Michely, L. Vazquez, 
and R. Cuerno, submitted. 

T. Hwa and M. Kardar, Phys. Rev. Lett. 62, 1813 (1989); 
Phys. Rev. A 45, 7002 (1992). 

L. Ponson, D. Bonamy, and E. Bouchaud, Phys. Rev. 
Lett. 96, 035506 (2006). 

E. Bouchbinder, I. Procaccia, and S. Sela, Phys. Rev 
Lett. 95, 255503 (2005). 

E. Bouchbinder, I. Procaccia, and S. Sela, J. Stat. Phys. 
125, 1029 (2006). 

A. Keller, R. Cuerno, S. Facsko, and W. Moller, Phys. 
Rev. B 79, 115437 (2009). 

A. Hansen, J. Schmittbuhl, and G. G. Batrouni, Phys. 
Rev. E 63, 062102 (2001). 

R. Pastor-Satorras and D. H. Rothman, J. Stat. Phys. 
93, 477 (1998). 

M. Nicoli, M. Castro, and R. Cuerno, J. Stat. Mech.: 
Theor. Exp., P02036 (2009). 

M. Nicoli, R. Cuerno, and M. Castro, Phys. Rev. Lett. 
102, 256102 (2009). 

M. Nicoli, R. Cuerno, and M. Castro, J. Stat. Mech.: 
Theor. Exp., P10030 (2011). 

M. Nicoli, M. Castro, and R. Cuerno, Phys. Rev. E 78, 
021601 (2008). 

Y. Saito, Statistical Physics of Crystal Growth (World 
Scientific, Singapore, 1996). 

E. Moro and R. Cuerno, Phys. Rev. E 63, 036104 (2001). 
M. Schwartz and E. Katzav, J. Stat. Mech.: Theor. Exp., 
P04023 (2008). 

M. Henkel and M. Pleimling, N on- Equilibrium Phase 

Transitions: Ageing and Dynamical Scaling Far from 

Equilibrium, Vol. 2 (Springer, Dordrecht, 2010). 

R. A. da Silveira and M. Kardar, Phys. Rev. E 68, 046108 

(2003). 

E. Katzav, Phys. Rev. E 68, 046113 (2003). 

E. Katzav and M. Schwartz, Europhys. Lett. 95, 66003 

(2011). 



